Cavity cooling of an optically levitated nanoparticle 
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The ability to trap and to manipulate individual 
atoms is at the heart of current implementations 
of quantum simulations [ , ], quantum comput- 
ing [^, ^], and long-distance quantum communica- 
tion [5—8]. Controlling the motion of larger par- 
ticles opens up yet new avenues for quantum sci- 
ence, both for the study of fundamental quantum 
phenomena in the context of matter wave inter- 
ference [ , ], and for new sensing and trans- 
duction applications in the context of quantum 
optomechanics [11, 12]. Specifically, it has been 
suggested that cavity cooling of a single nanopar- 
ticle in high vacuum allows for the generation of 
quantum states of motion in a room-temperature 
environment [ .^—^J\ as well as for unprecedented 
force sensitivity [16, 17]. Here, we take the first 
steps into this regime. We demonstrate cavity 
cooling of an optically levitated nanoparticle con- 
sisting of approximately 10^ atoms. The parti- 
cle is trapped at modest vacuum levels of a few 
millibar in the standing-wave field of an optical 
cavity and is cooled through coherent scatter- 
ing into the modes of the same cavity [^ ^, ^*^]. 
We estimate that our cooling rates are sufficient 
for ground-state cooling, provided that optical 
trapping at a vacuum level of 10~^ millibar can 
be realized in the future, e.g., by employing ad- 
ditional active- feedback schemes to stabilize the 
optical trap in three dimensions [20—2']. This 
paves the way for a new light-matter interface en- 
abling room-temperature quantum experiments 
with mesoscopic mechanical systems. 

Cooling and coherent control of single atoms inside 
an optical cavity are well-established techniques within 
atomic quantum optics [94-9^]. The main idea of cav- 
ity cooling relies on the fact that the presence of an 
optical cavity can resonantly enhance scattering pro- 
cesses of laser light that deplete the kinetic energy of the 
atom, specifically those processes where a photon that is 
scattered from the atom is Doppler-shifted to a higher 
frequency. It was realized early on that such cavity- 
enhanced scattering processes can be used to achieve 
laser cooling even of objects without exploitable inter- 
nal level structure such as molecules and nanoparticles 
[18, 19, 29, 30]. For nanoscale objects, cavity cool- 
ing has been demonstrated in a series of recent exper- 
iments with nanobeams [31-33] and membranes of nm- 
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scale thickness (e.g. [^/^^ '^^]). To guarantee long inter- 
action times with the cavity field these objects were me- 
chanically clamped, which however introduces additional 
dissipation and heating through the mechanical support 
structure. As one consequence, quantum signatures have 
thus far only been observed in a cryogenic environment 
[36, 37]. Freely suspended particles can circumvent this 
limitation and allow for far better decoupling of the meso- 
scopic object from the environment. This has been suc- 
cessfully implemented for atoms driven at optical fre- 
quencies far detuned from the atomic resonances, both 
for the case of optically trapped single atoms [26, ] and 
for clouds of up to 10^ ultracold atoms [38-40]. In con- 
trast to such clouds, massive solid objects provide access 
to a new parameter regime: on the one hand, the rigidity 
of the object allows to manipulate the center-of-mass mo- 
tion of the whole system, thus enabling macroscopically 
distinct superposition states [14, 15, 41]; on the other 
hand the large mass density of solids concentrates many 
atoms in a small volume of space, which provides new 
perspectives for force sensing [ 6, 17]. In our work, we 
have now extended the scheme to dielectric nanoparticles 
comprising up to 10^ atoms. By using a high-finesse op- 
tical cavity for both optical trapping and manipulation 
we demonstrate, for the first time, cavity-optomechanical 
control, including cooling, of the center-of-mass (CM) 
motion of a levitated solid object without internal level 
structure. 

To understand the principle of our approach, con- 
sider a dielectric spherical particle of radius r smaller 
than the optical wavelength A. Its finite polarizability 

^ = AneoT^Re < ^^ >(e: dielectric constant; eo: vacuum 

permitivity) results in an optical gradient force that al- 
lows to trap particles in the intensity maximum of an op- 
tical field [ ] . The spatial modes of an optical cavity pro- 
vide a standing- wave intensity distribution along the cav- 
ity axis X. A nanoparticle that enters the cavity will be 
pulled towards one of the intensity maxima, located a dis- 
tance xo from the cavity center. For the case of a Gaus- 
sian (TEMOO) cavity mode, the spatial profile will result 
in radial trapping around the cavity axis, hence provid- 
ing a full 3D particle confinement. In addition, Rayleigh 
scattering off the particle into the cavity mode induces 
a dispersive change in optical path length and shifts the 

cavity resonance frequency by Uo(xo) = ^^7/^ (l + ^) 

[ ]{^cav'' cavity frequency; Vcav'- cavity mode volume; 
xr: cavity- mode Rayleigh length). This provides the 
underlying optomechanical coupling mechanism between 
the CM motion of a particle moving along the cavity axis 
and the photons of a Gaussian cavity mode. The result- 
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FIG. 1: Optical trapping and readout of a nanoparticle in a Fabry-Perot cavity, (a) Nanoparticle in a cavity. A 

photo of our near-confocal Fabry-Perot optical cavity (OC) (F=76000; L = 2F%r =10-97 mm, determined via the free spectral 
range FSR). The white-shaded areas indicate the curvature of the cavity mirrors. The optical field between the mirrors traps 
a nanoparticle. The enlarged inset shows light scattered by the nanoparticle. (b) Schematics of two-mode optical trap 
and dispersive coupling. Two optical fields form standing-wave intensity distributions along the optical cavity axis (dashed 
lines; blue: control beam; red: trapping beam). Because of their different frequencies, the intensity maxima of the two fields are 
displaced with respect to each other. A nanoparticle is trapped at the maximum of the total intensity distribution (purple solid 
line). Since the trapping beam is more intense than the control beam, the nanoparticle is trapped at a distance x 7^ away from 
the control-beam intensity maximum xq. As a consequence, the nanoparticle oscillates within a region where the control-beam 
intensity varies with the particle position (blue arrow), resulting in linear dispersive coupling (see main text and appendix). 
The displacement x depends on the ratio between the intensity maxima of the two fields (c) Experimental setup. A Nd:YAG 
laser (A = 1064 nm) is split into three beams at the polarizing beam splitters PBSl and PBS2 (for simplicity waveplates not 
shown in the figure). The transmitted beam is used to lock the laser to the TEMOO mode of the OC and provides the trapping 
field for the nanoparticle. The beam reflected at PBSl is used to prepare the control beam, which is frequency-shifted by 
6uu close to the adjacent cavity resonance of the TEMOO mode, i.e., ^cj^FSR + A (A: detuning from cavity resonance). The 
single-frequency sideband at 6(jJ is created using an electro-optical modulator (EOM) followed by optical amplification in fiber 
(not shown) and transmission through a filtering cavity (FC) with an FWHM linewidth of 27rx500 MHz. The control and 
trapping beams are overlapped at PBS3 and transmitted through the OC with orthogonal polarizations. The OC is mounted 
inside a vacuum chamber (VAC). When a nanoparticle is trapped in the optical field in the cavity, its center-of-mass (CM) 
motion introduces a phase modulation on the control beam. To detect this signal, we perform interferometric phase readout of 
the control beam: At PBS4 the trapping beam is separated from the control beam and overlapped with the local oscillator (LO). 
After rotating the polarization, the control beam and the LO are mixed at PBS 5. High-frequency InCaAs photo detectors PDl 
and PD2 detect the light in both output ports of PBS5. We mix the difference signal of the two detectors with an electronic 
local oscillator of frequency FSR+ A and record the noise power spectrum of resulting signal at a spectrum analyzer (SA) (see 
Methods). 
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where we have allowed for a mean displacement x of the 
nanoparticle with respect to the intensity maximum xq 
{x\ CM position operator of the trapped nanoparticle; 
k=^: wavenumber of the cavity light field; h: cavity 
photon number operator). For the case of a single optical 
cavity mode, the particle is trapped at an intensity maxi- 
mum {x = 0) and, for small displacements, only coupling 
terms that are quadratic in x are relevant \3^. Linear 
coupling provides intrinsically larger coupling rates and 
can be exploited for various quantum control protocols 
[44]. However, it requires to position the particle outside 
the intensity maximum of the field. This can be achieved 
for example by an optical tweezer external to the cavity 
[14], by harnessing gravity in a vertically mounted cavity 
[45] or by using a second cavity mode with longitudinally 
shifted intensity maxima [13, 14]. 

We follow the latter approach and operate the optical 



cavity with two longitudinal Gaussian modes of different 
frequency, namely, a strong "trapping field" to realize a 
well-localized optical trap at one of its intensity maxima, 
and a weaker "control field" that couples to the particle 
at a shifted position x ^ 0. For localization in the Lamb- 
Dicke regime (/c^(f^) <C 1) this yields [ , ] linear op- 
tomechanical coupling between the trapped particle and 

the control field at a rate go = /7o(xo) sin(2A:x)/cw ^^^ 

per photon (m: nanoparticle mass; Qq: frequency of CM 
motion). Detuning of the control field from the cavity 
resonance by a frequency A = uj^^^ — uJc{oJc • control 
field frequency) results in the well-known dynamics of 
cavity optomechanics [12]. Specifically, the position de- 
pendence of the gradient force will change the stiffness of 
the optical trap, shifting Qq to an effective frequency (^eff 
(optical spring), and the cavity-induced retardation of 
the force will introduce additional optomechanical (pos- 
itive or negative) damping on the particle motion. From 
a quantum-optics viewpoint, the oscillating nanoparti- 
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FIG. 2: Experimental characterization of the nanopar- 
ticle cavity trap . (a) Schematic of the trap configura- 
tion. An optical cavity of length L — 10.97 mm is driven on 
resonance of a Gaussian TEMOO cavity mode by a laser with 
a wavelength of A = 1064 nm. The nanoparticle is optically 
trapped at position xq. Its center-of-mass motion in the axial 
direction of the cavity is described by a harmonic oscillator 
with a frequency Qo and an amplitude of approximately 10 
nm . In addition, the nanoparticle experiences collisions with 
the surrounding gas resulting in a damping rate 70. (b) Me- 
chanical damping 70 as a function of pressure. The solid 
line is a fit of kinetic gas theory to the data (see appendix D). 
(c) Position-dependent trapping frequency. The waist 
of the optical mode expands from approximately 41 /am at 
the cavity center to 61/jm at the cavity mirrors, resulting in 
a position-dependent trapping potential. Here, we show the 
corrsponding change of the trapping frequency ^^o with the 
position of the nanoparticle. (d) Power-dependent trap- 
ping frequency. We experimentally show the dependence 
of the trapping frequency on the intracavity power Pt. The 
solid lines in Fig. c, d are based on the theoretical model as 
described in the main text, with a scaling factor as the only 
free fit parameter. 



cle scatters photons into optical sidebands of frequencies 



ujc ± ^0 at rates A± 



9o{n)i^ 



known as Stokes 



4 («:/2)2 + (A±Qo 

and anti-Stokes scattering, respectively {n: FWHM cav- 
ity linewidth). For A > (red detuning) anti-Stokes 
scattering becomes resonantly enhanced by the cavity, 
effectively depleting the kinetic energy of the nanoparti- 
cle motion via a net laser-cooling rate of F = A- — Aj^. 
In the following, we demonstrate all these effects experi- 
mentally with an optically trapped silica nanoparticle. 

As is shown in Figure 1, our setup comprises a high- 
finesse Fabry- Perot cavity (Finesse F = 76000; n = 
2tt X 180 kHz) that is mounted inside a vacuum cham- 
ber kept at a pressure between 1 and 5 mbar. Airborne 
silica nanoparticles (specified with radius r = 127 ± 13 
nm) are emitted from an isopropanol solution via an ul- 



trasonic nebulizer and are trapped inside the cavity in 
the standing wave of the trapping field (see Methods 
Section). To achieve the desired displacement between 
the intensity maxima of trapping field and control field 
(x ^ 0), we use the adjacent longitudinal cavity mode for 
the control beam, i.e. the cavity mode shifted by approx- 
imately one free spectral range FSR =21;^ 13.67 GHz 
in frequency from the trapping beam (c: vacuum speed of 
light; L: cavity length). Depending on the distance from 
the cavity center xq, the two standing- wave intensity dis- 
tributions are then shifted with respect to each other by 
^(xo + L/2) (Figure Ic). For example, to achieve max- 
imal coupling qq for weak control beam powers, i.e. for 
/i = ^ <C 1 (^c(t)* Power of control (trapping) beam 
in the cavity), the nanoparticle needs to be positioned 
at xq = L/4, where the antinodes of the two beams are 
separated by A/8 [13, 14]. Note that when the control 
beam is strong enough to significantly contribute to the 
optical trap (/i ^ 0.1), the displacement x and both Qq 
and go are modified when /i is changed [39] . The ex- 
act dependence of these optomechanical parameters on /i 
depends on xq (see appendix A and [47, 48]). 

The optomechanical coupling between the control field 
and the particle can be used to both manipulate and de- 
tect the particle motion. Specifically, the axial motion 
of the nanoparticle generates a phase modulation of the 
control field, which we detect by heterodyne detection 
(see methods section). We reconstruct the noise power 
spectrum (NFS) of the mechanical motion by taking into 
account the significant filtering effects exhibited by the 
cavity (arising from the fact that k, ^ Qq) on the trans- 
mitted control beam ([49] and appendix A). The inferred 
position sensitvity of our readout scheme for a nanopar- 
ticle of approx. 170 nm radius is 4 pm/V^Hz, which is 
likely limited by classical laser noise (see below). 

The properties of our optical trap are summarized in 
Figure 2. The influence of the control beam on the 
trapping potential is purposely kept small by choosing 
/i ^ 0.1 and A ?^ 0. We expect that the axial mechanical 
frequency Qq depends both on the power of the trapping 
beam Pt and on xq through the cavity beam waist W{xo) 

via Qq = 
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ment with our data. The damping 70 of the mechanical 
resonator is dominated by the ambient pressure of the 
background gas down to a few millibar (Fig. Ic). Be- 
low these pressures the nanoparticle is not stably trapped 
anymore, while trapping times up to several hours can be 
achieved at a pressure of a few millibar. This is a known, 
yet unexplained phenomenon [21, 22, 47]. Reproducible 
optical trapping at lower pressure values has thus far only 
been reported using feedback cooling in three dimensions 
for the case of nanoparticles [21, 22] or, without feedback 
cooling, with particles of at least 20 /im radius [ ] . 

We finally demonstrate cavity-optomechanical control 
of our levitated nanoparticle. All measurements have 
been performed with the same particle for an intra-cavity 
trapping beam power Pt of approx. 55 W and at a pres- 
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FIG. 3: Cavity-optomechanical control and cooling of a nanoparticle. We obtain noise power spectra (NFS, see panel 
(a)) of the nanoparticle 's center-of-mass motion for different settings of the control-beam power Pc and detuning A. During 
each measurement, fi — ^ was kept constant ( Pt'. trapping beam power). Based on these NFS, we determine the effective 
mechanical frequency f^eff and linewidth 7eff of the optomechanical system, and its effective temperature Teff. We study the 
modification of these spectra caused by optomechanical interaction in panels (b), (c) and (e). Based on the data in panel (b) we 
infer the power- dependent strength of optomechanical coupling in panel (d). (a) Mechanical noise power spectra. Shown 
are examples of the mechanical NFS measured for constant control-beam power (/i = 0.3) at three different detunings A with 
respect to the cavity resonance frequency. The detuning results in a significant modification of the NFS due to optomechanical 
effects. Note that scale is changed by a factor of 5 in the bottom plot in panel (a). In order to determine the effective mechanical 
frequency f^eff and linewidth 7eff of the optomechanical system, we fit the NFS of an harmonic oscillator (black solid lines) 
to this data. We infer the value of the effective temperature Tgff from the equipartition theorem via direct integration of the 
NFS (see appendix C). (b) Optical spring. When the control beam is red-detuned from the cavity resonance (A > 0), we 
observe a characteristic modification of the mechanical frequency Qeff- The solid lines in (b) correspond to a theoretical model 
that is fitted to the data for each value of jjl. The optomechanical coupling go\/{fi) is one of the fit parameters (see appendix 
C). Based on these results for the optical spring, we calculate the theoretical expectations for 7eff and Teff, which are shown as 
dashed lines in panels (c) and (e). (c) Optomechanical damping. Linewidth broadening of the mechanical resonance as a 
function of the detuning A. (d) Optomechanical coupling. We infer the optomechanical coupling rate goy^Jn) from the 
strength of the optical spring (panel (b)) and show its dependence on the power ratio /i. This relation depends on the position 
xq of the nanoparticle in the cavity. For the data presented here, we determine xq = 1.56 ± 0.14 mm (see appendix E). We 
find very good agreement between the data and the theoretical model, where only the nanoparticle polarizability serves as a 
fit parameter (solid line; also see appendix C). (e) Cavity cooling. The decrease in effective temperature Teff is shown for 
increasing control-beam power. To obtain a good estimate of the measurement error, we average over measurements taken for 
detunings between A = 100 — 150 kHz (see appendix C). The dashed line is a theoretical prediction based on the parameters 
obtained from the fit to the optical spring data (panel (b)). 



sure of p ~ 4 mbar. This corresponds to a bare me- 
chanical frequency l]o/27r = 165 ±3 kHz and an intrinsic 
mechanical damping rate 7o/27r = 7.2 ± 0.8 kHz, respec- 
tively. Figure 3a shows the dependence of a typical noise 
power spectrum (NFS) of the particle's motion upon de- 
tuning of the control field. Note that the power ratio fi 
between trapping beam and control beam is kept con- 
stant, which is achieved by adjusting the control-beam 
power for different detunings. The amplitude scale, as 
well as the temperature scale in Figure 3e, is calibrated 



through the NFS measurement performed close to zero 
detuning (A = 1 kHz; blue NFS in Fig. 3a by using 
the equipartition theorem for T = 293K. This is justi- 
fied by an independent measurement that verifies ther- 
malization of the center of mass (CM) mode at zero de- 
tuning for our parameter regime (see appendix D). Both 
the inferred effective mechanical frequency r^eff (Figure 
3b and the effective mechanical damping 7eff (Figure 
3c show a systematic dependence on the detuning A 
of the control beam, in good agreement with the ex- 



pected dynamical backaction effects for linear optome- 
chanical coupling (see appendix A). A fit of the expected 
theory curve to the optical spring data allows estimat- 
ing the strength of the optomechanical coupling for dif- 
ferent values of /i (Figure 3d). If the position xq of 
the nanoparticle in the cavity is known, then this be- 
haviour is uniquely determined by Uo{xo). For a par- 
ticle position xq = 1.56 ± 0.14 mm, which was de- 
termined independently with a CCD camera, we find 
^o(^o) = 27r X (145 ± 2) kHz. These values allow to 
infer a nanoparticle displacement x ~ 0.15 x (A/2) = 77 
nm, yielding a fundamental single-photon coupling rate 
go ^ 27r X 1.2 Hz (for /i ^ 0). Assuming a (supplier 
specified) material density of p = 1950 g/cm and a di- 
electric constant esiOs =2.1, our results indicate a single 
trapped nanoparticle of radius r ~ 169 nm. 

The red-detuned driving of the cavity by the control 
laser also cools the CM motion of the levitated nanopar- 
ticle through coherent scattering into the cavity modes. 
Figure 3e shows the resulting effective temperature as 
deduced from the area of the NFS of the mechanical 
motion by applying the equipartition theorem. The ex- 
perimental data is well in agreement with the expected 
theory for cavity cooling (see appendix A). We achieve 
cooling rates of up to F = 27r x 49 kHz a nd ef fective op- 
tomechanical coupling rates of up to go \/ {fic) = 27r x 66 
kHz ((nc): mean photon number in control field), com- 
parable to state-of-the-art clamped mechanical systems 
in that frequency range [12]. The demonstrated cool- 
ing performance, with a minimal CM-mode temperature 
of 64 ± 5 K, is only limited by damping through resid- 
ual gas pressure that results in a mechanical quality of 
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25. Recent experiments [21, 2!] impressively 



demonstrate, that lower pressures can be achieved when 
cooling is applied in all three spatial dimensions. Given 
the fact that our cavity-induced longitudinal cooling rate 
is comparable to the feedback cooling rates achieved in 
those experiments, a combined scheme should eventually 
be capable of performing quantum experiments at mod- 
erately high vacuum levels. For example, our cooling rate 
is in principle sufficient to obtain cooling to the quantum 
ground state of the CM-motion starting from room tem- 
perature with a longitudinal mechanical quality factor 
of Q ^ 10^, i.e., a vacuum level of 10~^ mbar. Such a 
performance is currently out of reach for other existing 
cavity optomechanical systems with comparable frequen- 
cies. In addition, even larger cooling rates are expected 
when both beams are red-detuned to cooperatively cool 
the nanoparticle motion [ ]. 

Our experiment constitutes a first proof of concept 
demonstration in that direction. We envision that once 
this level of performance is achieved levitated nanopar- 
ticles in optical cavities will provide a room-temperature 
quantum interface between light and matter, along the 
lines proposed in [jlo, x^t, -t^:, ], with new opportuni- 
ties for macroscopic quantum experiments in a regime 
of large mass [15, 41, 52]. The large degree of optome- 
chanical control over levitated objects may also enable 



applications in other areas of physics such as for preci- 
sion force sensing [l ^^, l ] or for studying non-equilibrium 
dynamics in classical and quantum many-body systems 
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Methods 
Loading of nanoparticles into the optical cavity trap 

For our experiment we use silica nanospheres (Corpus- 
cular Inc.) with a radius of r = 127 ± 13 nm , which 
are provided in an aqueous solution with a mass concen- 
tration of 10%. We dilute the solution with isopropanol 
to a mass concentration of 10~^ and keep it for approx- 
imately 30 min in an ultrasonic bath before usage. To 
obtain airborne nanoparticles, an ultrasonic medical neb- 
ulizer (Omron Micro Air) emits droplets from the solu- 
tion with approximately 3/im size [ , ]. On average, 
the number of nanospheres per droplet is then approxi- 
mately 5 • 10~^. 

The nanospheres are loaded into the vacuum cham- 
ber by spraying the droplets through an inlet valve at 
the end of a 6mm thick, 90cm long steel tube. We keep 
the pressure inside the vacuum chamber between 1 and 
5 mBar via manual control of both the inlet valve con- 
nected to the nebulizer and the outlet valve connected 
to the vacuum pumps. During the loading process, the 
trapping laser is kept resonant with the cavity at the 
desired intracavity power for optical trapping. The low 



pressure minimizes pressure-induced fluctuations of the 
optical path length, which signiflcantly simplifles locking 
the laser to the cavity. 

Trapping in the conservative potential of the standing- 
wave trap is only possible with an additional dissipative 
process, which is provided fully by damping due to the 
remaining background gas. Within a few seconds af- 
ter opening the valve, nanospheres get optically trapped. 
The standing-wave conflguration provides multiple trap- 
ping positions. Trapped nanoparticles are detected by a 
CCD-camera, which is also used to determine their po- 
sition xq (see appendix E). If initially more than one 
position in the cavity is occupied, blocking the trapping 
beam for short intervals allows loosing surplus particles 
for our measurements. To move the trapped particle to 
different positions along the cavity, we blue-detune the 
control laser to heat the CM degree of freedom of the 
particle. The "hot" particle moves across the standing 
wave until the control beam is switched off and the par- 
ticle stays trapped at its new position (see Figure 2b). 



Readout of control beam 

For the position readout of the nanoparticle motion, 
we rely on the dispersive interaction with the control-fleld 
cavity mode. The control laser beam is initially prepared 
with a frequency difference of (^cj ~ 27r x 13.67 GHz with 
respect to the original laser frequency uq. When the 



control beam is transmitted through the cavity, it expe- 
riences a phase shift according to its detuning from the 
resonance oJcav Because the particle position in the cav- 
ity modifles the cavity resonance frequency uJcav^ a phase 
readout of the transmitted control beam allows recon- 
structing the nanoparticle's motion. To detect the phase 
modulation introduced by the particle motion along the 
cavity, we flrst mix the control beam with a local oscil- 
lator (LO, 3.15 mW; control beam power < O.lmW) at 
frequency uq at PBS5 (Fig. 1). In the output ports of 
PBS5, we then detect the optical signal at photodetec- 
tors PDl and PD2 (Discovery Semiconductor Inc. DSC- 
R410), which are fast enough to process the beat sig- 
nal at frequency Su. Their difference signal /(t), i.e., 
the heterodyne measurement outcome, contains the beat 
signal, whose phase 0opt is determined by the unknown 
path difference between the LO and the control beam. 
The beat signal carries sidebands representing the am- 
plitude and phase modulation imprinted on the control 
beam by the optomechanical system. We demodulate l{t) 
with an electronic local oscillator (ELO) with frequency 
Suj and phase ^elo (relative to the beat signal). From 
the resulting signal Sopt(^), we extract the phase modu- 
lation of l{t) by adjusting ^elo such that the total phase 
0ELO + ^opt = 7r/2 . This is achieved by locking the DC 
part of (sopt(^)) to zero. We record the NPS of 5opt(^) 
with a spectrum analyzer, which allows reconstructing 
the NPS of the nanoparticle's motion in post processing. 



Appendix A: Description of the Optomechanical System 
Optomechanical Hamiltonian 



To describe our experiment theoretically, we consider a nanoparticle that is optically trapped within a Fabry-Perot 
cavity. Two laser beams drive adjacent TEMqo cavity modes. One beam is used for optical trapping (trapping beam), 
the other for optomechanical control and readout of the nanoparticle center-of-mass motion (control beam). The two 
mode's resonance frequencies differ by one FSR = ^ {L: cavity length). In the most general case, the two lasers can 
be detuned from the respective cavity resonance frequency by At and Ac {^c{t) • detuning of the control (trapping) 
beam). The system is described using the following Hamiltonian [4b]: 



H/h = AtdlcLt + Acajttc 



2mh 



;it 



Uodlat sin {ktx) 



-Uodldc sin^(/ccx) + iEt{d\ — dt) + iEc{d\ — dc)-, 



(Al) 



where Uq can be understood as the cavity resonance frequency shift introduced by a nanoparticle that is located 
at the intensity maximum at the center of the optical cavity. At the same time, HUq is also the trap depth created 
by a single intracavity photon {dK^. / dc^t)' creation/annihilation operator of the control (trapping) fleld in the 
cavity; m: mass of the nanoparticle; x (pm)'- position (momentum) operatior of the nanoparticle's CM; kc{t^/E(,(j^y. 
wavenumber/driving fleld of the control (trapping) beam). 

Given |/ct — /cd <C /Cc, one can regard {kt — kc)x as a position-dependent phase shift between the standing waves of 
the two intracavity flelds via sin^ (ktx) = sin^ {kcX -\- {kt — kc)x) = sin^ {kcX -\- (/?), where 



cp = {kt- kc)Q 



27rFSR 



We include the dependence on kt in cp and use k instead of kc from this point on. Further, we rewrite the position 
operator x as the sum of three terms: x = Xq -\- x -\- f ^, where Xq is the position of the intensity maximum of the 



control field with respect to the cavity mirror (x: the nanoparticle's mean displacement from Xq, Xm- the nanoparticle's 
displaced position operator with (xm) — 0). Note that in the main text we always use the distance from the cavity 
center xq, where Xq = xq + L/2. We also introduce the dimensionless position operator ^x with x^ — -^gs • ^^7 where 
bx = -j^(b-\-b^) {Xgs'- Ground state extension of the mechanical oscillator, b^^"^: CM-motion annihilation (creation) 
operator). 

We approximate the trigonometric functions in equation Al to a second-order in Xm and perform a displacement 
operation of the light operators: aj -^ aj + clj about their steady-state mean values at and ac- The Hamiltonian 
after these modifications is: 



H 



A^l^tp + Ac|acp + Atat{at + aJ) + Acadctc + 4) + ^t^I^t + ^cCtlcLc 
1^ - Uo\atf sin2(fc(4 + x) + (^) - Uo\a,\^ sm\k{x'^ + x)) 



— 2Uok^\atf cos(2/c(4 ^x)^2ip)-^- 2Uok^\ac\^ cos(2/c(4 + x))-^ 

— Uokat sin{2k{xQ -\- x) -\- 2(p){at + ^l)^m — U^kac sin(2/c(xo + x)){ac + ^J)f ^ 

— /7o|at|^^sin(2/c(xo + x) + 2Lp)xra — /7o|ac|^^sin(2/c(xo + x))xra 

— Uoat sin^(A:(xQ -\- x) -\- (p){dt + aJ) — Uoac sin^(/c(xQ + x)){ac + ^J) 

-h iE'tlaJ - at) + i£;c(4 - dc). (A2) 

Line A2 takes the form of a harmonic potential ^ ^^"^ with mechanical frequency Qq: 

2hU k'^ 

Ql = ^— {\at\^ cos{2k{xQ + x) + 2(p) + \acf cos{2k{xQ + x))) . (A3) 

Line A2 determines the linear dispersive coupling of the nanosphere CM motion to the trapping and cooling beam. 
Note that the trapping beam also shows linear coupling when the cooling beam is strong enough to significantly 
contribute to the optical trap: 

go,t = Ct • ^gs = UokXgs sin 2{k{xQ + x) + (/:?) 

go,c = Cc • ^gs = UokXgs sm2k{xQ + x). (A4) 

Note that we use go = ^o,c in the main text. To study the dynamics of the system, we solve the Langevin equations 
for both light fields: 

dt = -{- ^i{At-Uosm'^{k{xQ^x) ^(f)){dt^at) ^ Et - CtO^tXm 

dc = -(- +i(Ac - /7osin^ A:(xo + x))(ac + ac) + £^c - CcCi^c^m (A5) 

The additional loss terms account for the cavity amplitude decay rate hi. The value of k, is assumed to be equal for 
both light fields due to the small difference in their wavelengths. For the steady-state solutions of dt and dc we find: 

Et 

at = 



+ i{At — Uq sin k{xQ -\- x) -\- (p) 



2 



f + i (A, - Uo sm\k{x'Q + x))) ' ^^^^ 

In our experiment, the Pound-Drever-Hall feedback loop keeps the trapping-laser frequency resonant to the corre- 
sponding cavity resonance frequency when the particle is in its steady state position. In other words, the detuning 
At compensates the frequency shift caused by the particle such that A^ — Uq sin^(/c(xo -\- x) -\- (p) =0. 

On the other hand, the frequency of the control beam is varied throughout the experiment. We are interested in the 
detuning A of the control beam with respect to the cavity resonance when the nanoparticle is located at its steady 
state position: A = Ac — Uq sin^ k{xQ -\- x). 

The trapping beam power is not changed throughout the experiment. In contrast, the control beam power is always 
set to achieve the desired ratio between the power of the two intracavity fields /i: 

., _ |a.|^ _ m' (f )^ 

12 ~ 117 12 /„x2 . „■ V^V 



K|2 |S*|2(|)^ + A2 



Heisenberg's equation of motion for the particle becomes: 

Xm ^ftl{fi)xm = \\at\'^ sin 2{k{xQ + x) + (/?) + \ac\'^ sin 2k {xq + x)1 - 7m^m (A8) 

m 

where we included an additional damping term 7m^m, which is due to the collisions of the nanoparticle with the 
surrounding gas. 

From Equation A8 we find a steady state condition on Xq + x, that enables us to determine the mechanical frequency 
fto and the displacement x as a function of /i: 



nlifi) = ^o(0)\/l + M^ + 2/icos2(/:) (A9) 

_ _ sin2(f 
tanz/cx = . 

/i + cos 2(p 

Thereby, the mechanical frequency in absence of the cooling beam is (equation A3) : 

. 2hUoe , 2 

^o(0)= ^ \(^t\ . 

Note that the case of a control beam that significantly contributes to the optical trap that has been presented here 
has also already been published in [47, 48]. 

Cavity mode shape 

Up to this point, we have neglected the mode shape of the TEMOO cavity mode (Fig. 2a, main text). The waist of 
the mode, however, depends on the position xq in the cavity. The maximum intensity of the standing wave along the 
TEMOO mode in the cavity is, accordingly, position dependent [55]: 

I{xo) = /o -^, 



note, that we have used here xq as the distance from the center of the cavity. It is related to the distance from the 
mirror Xq hy xq = Xq — ^ {xr: Rayleigh length of the mode). Therefore, Uq is an explicit function of the trap position 

Xq: 

{oJcav'- laser frequency, ^: particle polarizability, Sq: vacuum permittivity, Vc'. cavity mode volume). The polarizability 
of a particle is (see e.g. [x^]): 

/e — 1 



{e: nanoparticle's dielectric constant; r: particle radius). In the main text we use these equations to determine the 
estimated particle size from Uo{xo), which is determined from the control beam power dependent coupling ^o(see 
main text, figure 3d) and the independently determied position of a particle in the cavity xo(see appendix E). 

Langevin Equations, effective frequency and damping 

The Langevin equations for the mechanical quantum harmonic oscillator coupled to a thermal bath are: 

A _ Pm 

m 

Pm = -mftlXm - ImPm + ^ hCjaj{d] + ttj) + 7^(t), (AlO) 



where r]{t) is a thermal noise term, with the fohowing correlation property [56]: 



<■*'*'» =i^/l--'-'''--'>G£r- 



We assume that we are in a temperature range where ksT /h ^ VLq: 

For the light beams, we can use the equations of motion as provided in equation A5 after the displacement of the 
light operators: 

^c = -{'^^iA)dc^iCcacXm^\^{cf^df), (All) 

dt = - 2^^ + KtO^tXm + \/^{^r + ^t"") 

By Fouriertransformation, we obtain a linear system of equations from which we retrieve the final expression for 
the position spectrum Sxx{^) of levitating nanoparticles CM motion: 

^xx = \Xm I [^th + ^rpj^ 

where 5'th is the thermal noise contribution and ^rp is the radiation-pressure contribution. In the regime our experi- 
ment is currently operating (T = 293 K; air pressure approx. 1-5 mbar), we expect that the thermal-noise contribution 
pre vales: 



5'th = ^os7^ 



2 2kBT 



The effective susceptibility of the mechanical oscillator is 



(l]eff(^)^ - ^^)^ - i7eff(^)^^^ 






where, following [''^], we used the expressions: 

7eff(a;) = 7m - 7 5 ^7 5 V 

((f)V(c. + Ap)((f)V(c.-A)2) 



^eff(^) = 



v2 






.^2. 



((!)' + ('^ + A)2)((t)' + (c.-A)2) 



Appendix B: Position readout by homodyne detection of the control beam 

The expressions for the mechanical oscillator's dynamics, as well as its relationship to the control beam in the cavity, 
have been derived in the previous section (equations AlO and All). In the following two sections, we will discuss 
how the mechanical oscillator position NFS is determined from the NFS obtained by homodyning of the control-beam 
phase signal in transmission of the cavity (see Methods M2 for implementation of homodyne detection). 

We first derive the control light field in cavity transmission via the cavity input-output relation [57]: 

dl''\i) = ^dc(t) - d^^(t), (Bi) 

where d^^ describes the quantum noise at the cavity back mirror (i.e., the side from which the cavity is not driven). 
Even though our detection scheme occurs in two steps as described in the methods section, it is completely equivalent 
to a standard homodyne detection. The output signal is accordingly described by [38]: 

^opt(t) = \ {\dT' + dLo\^ - \dT' - dLo\^) = dT'dlo + dT'^dLo. (B2) 
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where we describe the local oscillator by aLo{t) = (^lo ' 6"**^^^^+^^ where determines the detected quadrature of 
the control beam and a^o is assumed to be real. In our experiment, the readout phase is locked to measure the phase 
quadrature: ^ = f • 

From equations Bl and B2 we obtain: 



(|5opt(^)l ) = ^Cc^c IXc(^) +Xc(-^)I Sxa:{^)S{uj), 



(B3) 



where we used the cavity susceptibility Xc(^) = ^_i(lj_/\) - Up to a proportionality factor, Equation B3 resembles the 
result of the detection described in [ ], and allows us to derive the mechanical NFS Sxx{^) from the detected signal. 

Appendix C: Data Evaluation and Temperature Calibration 

To extract the mechanical NFS, we first measure the spectrum of the homodyne phase readout with and without 
particle for all values of fi and A. We obtain (|5opt('^)P) by substracting the background NFS (without particle) 
from the NFS with particle. To reconstruct the mechanical NFS Sxx^ we need to account for the filtering by the 
Fabry-Ferot cavity. We therefore divide (|sopt(co')p) by \Xc{^) + Xc(~^)P following equation B3. The exact shape 
of Sxx is given by equation A 12. To determine the effective frequency, damping and temperature we assume that we 
can describe the CM-motion of the particle as an harmonic oscillator, which is fulfilled as we are not operating in the 
strong coupling regime: 



O Teff(from integration) 
• Teff (from fit) 

T theory 
eff 




100 150 200 250 300 
A [kHz] 



FIG. 4: CM-motion temperature as a function of detuning. The temperature of the CM-motion along the cavity axis as 
a function of detuning. The values are inferred via the equipartition theorem from the direct integration of the NPS (T/f^, solid 
circles) and from the fitted spectra (T*^, empty diamonds). The solid line shows the theoretical expectation T^^^°^'^(A) infered 
from the detuning dependent frequency fit (optical spring; main text, Fig. 3). 



f{uj) = a-Tef^' 



(Cj2 



^e^^- 



-^'7eff' 



(CI) 



By fitting this model to Sxx^ we obtain 7eff, ^eff and T^*^. The calibration constant a is determined such that 
Teff = 293K in a particular measurement that was performed close to zero detuning (A = 1 kHz for /i = 0.4, blue 
NPS in Fig. 3, main text). This results in the values for the optical spring l^eff and damping 7eff in Fig. 3, main text. 

We can determine the optomechanical coupling go from the the detuning dependence of l^eff for a given value of 
/i. However, we do not have an explicit analytical expression for this dependence. Instead, we apply the following 
strategy: 

Using equation A 12, we can calculate the optomechanical NPS S^^"^^^ of our system for a given set of parameters 
(/>:, ^0, ^0 and 5 A). Here, SA is a systematic deviation from the detuning we set in the measurement: each value of 
A can be set precisely up to the uncertainty in the actual cavity resonance frequency. This frequency difference is 
accounted for with a joint offset SA in the values of A that is used as a fit parameter. We treat 5'*^®°^^ in the same 
manner as the data and extract 7eff^°^"^, ^eff^°^"^ ^^^ ^eff^°^"^ ^y fitting f{uj) for each value of A. We use ft^f^'^^^{A) 
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FIG. 5: Measurement of the NPS-area as a function of the ambient pressure. The data were taken for a resonant control beam. 
The area is proportional to the temperature of the CM motion of the nanoparticle for a given value of /i. It is independent 
of pressure indicating that the CM motion of the trapped nanoparticles is thermalized with the surrounding gas over the full 
measurement range. We conclude that the effective temperature of the CM mode is room temperature (293 K). The scatter of 
the data corresponds to an standard deviation of 5% for a temperature that is infered from the area of the NFS. 



as a model that we fit to lleff optimizing the parameters go, Qq and Su in a least-square fit. The FWHM cavity 
line width tz is determined independently. The best fit parameters are used to obtain the theoretical dependences of 
%ff^°^"^ and ^eff^°^"^ ^^ ^^^ detuning shown in Fig. 3 in the main text and T^^^^^"^ shown in Fig. 4. 

The corresponding values of the predicted effective temperatures T^^^'^^^ are shown in Fig. 4 along with the ex- 
perimental data for /i = 0.4. The latter is obtained in two ways: firstly as a free parameter T^*^ in the fitted model 
f{uj) and secondly by direct integration over the measured NFS via T/^ = a^Vt^^ j Sxxd^- The calibration factor 
a^ is derived in the same way as a. The values T^*^, obtained via fitting, agree well with those obtained by direct 
integration of the NFS. For small detunings A, the data follows the theoretical curve, while for larger detunings, 
heating unaccounted for in the theoretical model seems to occur. We are still investigating this effect, which may 
be due to laser noise. To obtain a good estimate of the minimal temperature achieved experimentally, we average 
the temperature obtained for a range of detunings A/27r G [100, 150] kHz. The range is chosen such that the onset 
of temperature increase is not yet strong and the predicted range of Teff is small compared to the distribution of 
measured temperatures. The experimental data in Fig. 3e in the main text is obtained by applying this evaluation 
for the different values of ji for T/^ obtained by direct integration. The theory curve in Fig. 3e in the main text is 
obtained by averaging the theoretical prediction for T^^®°^^ over the same range of detunings A/27r G [100, 150] kHz. 



Appendix D: Kinetic gas theory - Pressure-dependent damping 

The pressure dependence of the damping for a trapped nanosphere is given by [21, 22, 59] : 

GTTT^r 0.619 



7o 



m 0.619 + irn 



(1 + Ck) 



(Dl) 



where 77 is the viscosity coefficient for air, r and m are the radius and mass of the nanosphere, Kn = Afp/r is the 
Knudsen number and Afp is the mean free path for air particles. Ck = 0.31i^n/(0.785 + 1.152Kn + Kn^) is a small 
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correction factor necessary at higher pressures [^^]. Figure 2b in the main text shows a pressure dependent damping 
measurement, where the control beam is used just for readout (i.e. /i=0.1, resonant). 

Figure 5 shows the temperature associated with the CM motion of the particle as a function of pressure. At high 
pressures, the nanoparticle experiences more collisions with the gas resulting in a stronger damping of its CM motion. 
If the nanoparticle CM motion was not thermalized at low pressures due to a heating process, better thermalization 
and therefore lower temperatures would be expected at higher pressures due to the increased damping rate. As Fig. 5 
shows a constant CM motion temperature for the different pressures and values of /i, we conclude it is thermalized 
with the environment in all these measurements, which implies a temperature of 293K as long as no optical damping 
is introduced. 




FIG. 6: Schematic of the configuration of the CCD imaging setup. We use a combination of three COD cameras to 
observe and locate the particle. This configuration allows access to a larger range of positions along the cavity axis compared 
to a single camera. Based on the pictures of the CCD cameras, we determine the position xq of a particle in the cavity. 



Appendix E: Position detection 



Three cameras with achromatic lenses monitor the cavity and image the light scattered off trapped nanoparticles. 
As can be seen in Figure la of the main text, the black retaining rings and the concave shape of the mirrors prevent 
optical access over the whole cavity length from a single point of view. We use a configuration of three CCD cameras, 
as shown in figure 6, to extend the field of view. By combining the images from the 3 cameras, we can reconstruct a 
larger field of view. To determine the position of the particle from the image, we need to calibrate the coordinates. To 
this end, the mechanical frequency at several positions is measured along with the position of a particle on the CCD 
image. This measurement is repeated for several particles. The frequency dependence on position allows calibrating 
the camera. 

The mode shape in the optical cavity is well-known from the curvature of the mirrors and the cavity length, which 
is determined with high precision from the FSR. The expected longitudinal frequency dependence of a nanoparticle 



trapped in the standing wave is Qq = Qc 



{xc'- cavity center position, Qc- frequency at position 



Xc] xr: Rayleigh length of the Gaussian mode). The measured mechanical frequencies for several different trap 
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positions of the same nanoparticle and the corresponding coordinates ( (in pixels) on the camera images are fitted 
to the function flo = Qc / "*" = with fit parameters (c (coordinate of the center of the cavity in pixels), 

^ (conversion factor between pixels and millimeters) and Qc (mechanical frequency in the center of the cavity). A 
corresponding measurement for one nanoparticle with calibrated length scale is shown in the main text, figure 2b. 
Based on this calibration we determine the position of the nanoparticle used in the measurements summarized in 
Fig. 3. It is located at a distance Xq = 3.92 ± 0.14 mm from the cavity mirror, i.e. at a distance xq = 1.56 ± 0.14 mm 
from the center of the cavity. 
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